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We give analytic upper bounds to the channel capacity C for transmission of classical informa- 
tion in electromagnetic channels (bosonic channels with thermal noise). In the practically relevant 
regimes of high noise and low transmissivity, by comparison with know lower bounds on C, our in- 
equalities determine the value of the capacity up to corrections which are irrelevant for all practical 
purposes. Examples of such channels are radio communication, infrared or visible-wavelength free 
space channels. We also provide bounds to active channels that include amplification. 



Shannon famously proved that the maximum num- 
ber of bits transmitted through a narrowband Gaussian- 
noise channel is C = log 2 (l + S/N) for each use of the 
channel, where S/N is the signal to noise ratio. At bot- 
tom, the noise has a quantum origin, and the calcula- 
tion of the capacity requires a quantum description of 
the channel. Accordingly, one of the oldest questions in 
quantum information theory is the calculation of channel 
capacities || . 
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FIG. 1: Plots of the bounds for the passive channel with 
thermal noise (a) and for the amplifier (b) in the practically 
relevant regimes of high thermal noise (large N): the upper 
and lower bounds (red and black curves of Figs. ^|b and " 
are basically coincident here, (a) Plots of the bounds ([!]), 
and of the bound of Koenig and Smith |Q| as a function of t : 
transmissivity r\ for a microwave radio communication, , 
with N = 2000 at room temperature, TV = 10 11 for ~ 1 mW 
transmission power (assuming that one channel use lasts for 
one oscillation period of the radiation). In this regime the 
lower bound (111) is indistinguishable from our upper bound 
(^) (black curve) (see also the magnification in the insert), 
whereas the gap with Koenig and Smith's bound (green curve) 
is evident. Typical transmissivities are very low for these 
channels (e.g., r\ = 0.04 for a 14 dB attenuation), (b) Capacity 
of the amplifier A% as a function of the gain k: here the 
lower bound (|l|) is indistinguishable from the upper bound 
(^) (black line), whereas the upper bound ( pi| ) (green line) is 
not useful. Here again N = 10 11 and N = 2000. 

Some of the most practically relevant communication 
channels are active and passive bosonic channels with 
thermal noise, e.g. radio or infrared-light communication. 



In this paper we provide upper and lower bounds for their 
capacity. In the case of the passive bosonic channel, our 
bounds supersede the recent one given by Koenig and 
Smith Q: in particular in contrast to their bound, for 
the practically relevant regimes of large thermal noise or 
low transmissivities (where each channel use can convey 
only small fractions of a bit) our bounds arc sufficiently 
tight to constitute an expression for the capacity which 
is good for practical purposes. These findings are consis- 
tent with the Holevo- Werner conjecture [||, that Gaus- 
sian mixtures of coherent states achieve capacity and that 
these channels are additive. In other words, in situations 
of practical interest, quantum effects (such as entangle- 
ment among subsequent channel uses) do not give any 
advantage and a coding alphabet composed of coherent 
states (e.g., the is output from a mascr or laser) achieves 
capacity. It is important to stress, however, that there 
are other regimes in which our inequalities are not tight: 
(slight) quantum advantages in low-noise regimes might 
be still possible. 

We consider two passive channels: the thermal bosonic 
channel £^ that can be modeled by a beam splitter of 
transmissivity 77 that mixes the signal with a thermal 
state with mean photon number N (the capacity of this 
channel for N = is already known H), and the classi- 
cal additive noise channel Af„ in which the signal is ran- 
domly displaced in the complex phase-space according 
to a Gaussian probability distribution of variance n. We 
also consider phase-insensitive amplifiers with gain 
K ^ 1, whose additional input mode (required to ensure 
the correct commutation relations of the fields) is in a 
thermal state of mean photon number N. Using a signal 
of N average photons, with an alphabet of Gaussian- 
distributed coherent states, one find the following lower 
bounds for their capacities C: 



C{£») > g( V N+(l- V )N)-g((l- V )N) , 

C(AQ > g(N + n)~g(n) , 

C(A%) > g( K N + ( K -l)(N + l))-g(( K -l)(N- 
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where g(x) 
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1) — xlog 2 x. (A proof 



of the Holevo- Werner conjecture would turn these into 
equalities, but it has been elusive even after a decade 
of concerted efforts §-§.) The main result of our pa- 
per is a collection of upper bounds which asymptotically 
match the above lower bounds in the practically relevant 
regimes of high noise, or low transmissivity, or high am- 
plification (see Fig. |l|). In particular we show that the 
following inequalities apply 

C{£») < 9 (vN + (1 - n)N) - ,g((l - V )N - v ) , (4) 
for T) < jfzj , 

C(M0 < g(N + n) - g(n - 1) for n ^ 1, (5) 
C(A%) < g(KN + ( K -l)(N+l))-g(N( K -l)-l) , 
for K>Z$±. (6) 

In addition to these simple bounds (which arc nonetheless 
good enough for many applications) we also derive other, 
even tighter, bounds in what follows. All our bounds 
apply to narrowband channels, but they can be extended 
to broadband channels using the variational techniques 
we detailed in [|w| . The remainder of the paper is devoted 
to the proof of these and of the further bounds. 

On Bounds and Conjectures: — To characterize the 
channels one can use their action on the state's char- 
acteristic function := Tr[p e^ a ~ M a ] (a being the 
annihilation operator of the mode): x(aO is transformed 
by the channels S^, Af n , and A% as (see, e.g., Ref. jfjj) 



x(m) 






(7) 
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x(Kfx) e -C«-i)C^+Va)lfl a . 





The main difficulty in the calculation of the classi- 
cal capacity of a quantum channel $ is superadditivity 
||: there exist channels in which the alphabet that 
achieves capacity must be composed by entangled quan- 
tum states that span multiple channel uses. Accordingly, 
one must regularize as follows || 



C($) 



lim C x (<f>® m )/m, 



(8) 



where indicates m uses of the channel 4>, and 12, 13[] 
C x (*) - max s(j2pMPi}) - £PiS(*[Pi]) ■ (9) 

Here S(p) := — Trplog 2 p is the von Neumann entropy, 
^[pi] is the output state from the channel VP (that may 
represent multiple uses of $), and the maximization is 
performed over the set of ensembles {pi,Pi} formed by 
density matrices pi and probabilities pi that may satisfy 
some resource constraint (such as on the average photon 
number N discussed above). Lower bounds to C($) can 



be obtained by calculating the right hand side of (||) for a 
specific encoding alphabet, i.e. fixing the value of m (say 
m = 1) and using a specific choice of for pi and pi, as 
was done to obtain the inequalities In contrast, 

an upper bound for C($) is provided by 



lim S„ 



l ($ 0m )/m 



(10) 



where S ma x(&) = max p S(<&(p)) is the maximum output 
entropy for a single channel use [using the same restric- 
tions in the maximization as in the definition of C($)], 
and Smin^) = min p S(^(p)) is the (unrestricted) mini- 
mum output entropy of the channel \&. The regulariza- 
tion over m in ( pij| ) is required by the superadditivity of 
the minimum output entropy [ |PT| , and constitutes the 
main difficulty in deriving bounds through (pX| ). How- 
ever, if $ = $bb is entanglement-breaking [p_4|, p"5[, the 
rcgularization is unnecessary |l(| and (|Io| ) can be re- 
placed by 

C{<$>EB) «S S max {$ EB ) - S m in(^EB) , (11) 

[notice that both S min ($%g) and C x ($|2) are additive 
quantities] . 

For i> = £^ , Af n , or A% the first term on the right of 
inequality (flO|), S max ($), is easily computed by exploit- 
ing the fact that the thermal state maximizes the entropy 
for fixed average photon number N: 



c (pis 

u max y-'-q 

S (A N 



g(r,N + (1 - rj)N) , 
g(N + n) , 

g( K N+( K -l)(N + l)) 



(12) 



In contrast, evaluating the second term on the right 
of inequality ( |Io[ ) is extremely demanding: the Holevo- 
Werner conjecture can be rephrased into a conjecture on 
the values of S min ($® m ) §, 0, |jj, which states that 
the min is achieved by a vacuum state |0)® m . If this 
were true, one could use ( |To| ) to provide upper bounds 
that exactly match the lower bounds (flj)-(||). A proof 
of this is lacking, but in Ref. || several bounds were 
obtained for the special case of rn = 1: they constrain 
S m in{£^ ) and S m i n (J\f n ) close to their conjectured values 
of g((l — rj)N) and g(n), respectively. Using (|ll|), such 
bounds can be immediately translated into constraints 
on the capacity C whenever the maps are entanglement- 
breaking, i.e. when n ^ N/(N + 1) for £!? , when n ^ 1 
for Nn, and when N ^ 1/(k - 1) for A" @. For in- 
stance, exploiting this fact, inequality (^ 
rived by replacing the term S m i n {£^) of 
single- mode lower bound A of Ref. p[ . 
ally, the same approach exploited in p| can be adapted 
to the multi-channel use scenario to construct tight in- 
equalities directly for the quantities 5 m i n ([£,^]® m )/m 
and S min {[Af n ]® m )/m. When substituted into @ to- 
gether with the identities (|l^) these then translate into 
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FIG. 2: Plots of the bounds in regimes that emphasize the gap 
between the lower and the upper bounds (these regimes are 
typically not interesting in practical applications), (a) Ca- 
pacity of the Gaussian channel M n for N = 1. Red curve: 
lower bou nd (H ); blue, black, green curves: upper bounds (|B|), 
([L6|), and (|l7|)7 respectively. The yellow area emphasizes the 
gap between the best upper and lower bounds, (b) Capac- 
ity of the passive electromagnetic channel £^ . Red curve: 
lower bound ([j]); blue curve: upper bound (Q) (valid only 
for r) ^ N/(N + 1) shown as a vertical dashed line where 
the channel becomes entanglement breaking); green curve: 
Koenig and Smith's bound from 0|; black line: upper bound 
(p7|); black dashed line: upper bound (p^). Here N = 0.5 
thermal photons and N = 1 average photons in the signal 
(which gives bits-per-photon for each channel use), (c) Plots 
of the bounds for the amplifying channel A^ with gain k. 
Red curve: lower bound (jij); black curve: upper bound ([]) 
(valid only for k ^ (TV + l)/iV); green curve: upper bound 
([||. The discontinuity for k = (N + 1)/N (vertical dashed 
line) separates the entanglement breaking regime on the right 
from the pure- loss regime on the left. Here N — 3 and N = 1. 



a collection of upper bounds for C that hold beyond the 
entanglement-breaking regime detailed above. 

Bounds for the Additive Classical noise channel Af n : — 
As detailed below, the bounds a, b, and d of Ref. || for 
m = 1 can be generalized to arbitrary m as follows 

S m ra{K m )lm > g(n - 1) , [Vn > 1] (13) 
S min {M® m )/m > log 2 (2n+l), (14) 
5 min (7V ! f ro )/m > l + Iog 2 (»), (15) 

whence, using @, Eq. Q gives (|), while Eqs. (0) 
and ( |l5| ) respectively give the further bounds 

C(AT n ) < 5 (Ar + n)-log 2 (2n + l), (16) 
C(AT n ) < ff (JV + n)-l-Iog 2 (n), (17) 

[the generalization of the bound c of B is not reported 
here since it converges to Eq. dTl) for m — > oo]. These 



bounds are compared to the lower bound (0) in Fig. ||(a): 
note how the gap between the upper and lower bounds 
closes asymptotically for high noise, n — > oo. 

The proof of Eq. ([ilf), and hence of the bound (||), was 
given in Ref. Jl7[ by expanding a generic input state p 
in terms of its multi-mode Husimi distribution function 
and applying the concavity of von Neumann entropy. An 
alternative proof follows from inequality a of Ref. and 
from (fTl|), using the fact that the channel Af n is entan- 
glement breaking for n > 1 [jl5] . 

The proof of Eq. ( |i"4| ) exploits the fact that the von 
Neumann entropy is never smaller than the Renyi en- 
tropy of order 2 |g |§ i.e. S(p) ^ S 2 {p) := - log 2 Tr[p 2 ]. 
Thus, for all input density matrices p of m channel uses 
we have 

W m (/>)) > SM m { 9 )) > m log 2 (2n + 1) , (18) 

where the last inequality follows from the fact that the 
minimum Renyi entropy of integer order at the output of 
the channel M n is additive and saturated by the vacuum 
input state The bound Q, and hence @, follow 
by minimizing with respect to p. 



The proof of Eq. ( 15 ) closely follows the proof of bound 
d in Ref. || for m = 1. Indeed, given a generic pure input 
state the eigenvalues "fk of the relevant output state 
p' = A/"® m (|'0)(V , |) can be expressed as 



7 ,= / d 2m pPi m \mik\D{m)t 



(19) 



where py^) is the corresponding eigenvector of p' , £>(/!) 
is the m-modc displacement operator, Pn m \p) := 
cxp[— |y(2| 2 /n]/(7rn) m , and the integral is performed over 
the m-dimensional complex vectors /7 G C m . By convex- 
ity, for all z 1 one can write 

d 2m r! 



Tr[(p'Y 



E 



z — l\m 



l/(zn z " 1 ) 



(20) 



which gives inequality (|15[), and hence (|l_7j), by remem- 
bering that S(p') = lim Ml+ log 2 Tr[(p') z ]/(i - z) E| 



Bounds for the Lossy Thermal channel £^ : 



The 

bounds (^3|)-(|15|) can be immediately turned into in- 
equalities for S m i n {[£^]® m ) by exploiting the compo- 
sitions rules Q that link and Af n that also ap- 
ply to the multi-use scenario m > 1. In particular, 



A^f™ o [£°]® m Hence, following the same 
, we find 

) ■ (21) 



reasoning of Ref. |6J , we find 



In particular, from ( |l4] ) and ( p"5| ) we obtain the multi-use 
versions of the bounds B and C of Ref. 0, i.e., 



-'mm 



([£^r n )/m > log 2 (2(l- ^jV + l) 
„([^]® m )m ^ l + log^l-^JV) , 



(22) 
(23) 
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which give rise to the following bounds for the capacity 

C(S») < g((l - + N) - log 2 (2(l - V )N + 1) , 

(24) 

C(£?) < g((l - r))N + N) — 1 — log 2 ((l - r,)N) , 

(25) 

[the bound obtained from (|l3| ) is not reported here as it 
is always subsided by the inequality (0)]. Further bounds 
can be obtained by generalizing the inequalities E and F 
of Rcf. ||: for all integers fc, we have for inequality E: 

*J min\Y 1 -' rj J / 



to 

nin ([C ie 



> ^ 9(^(1 -V)N) for V ^l/k, 



m 



(1-7^-77+- 



for r/ 1/fc, 



(26) 



and for inequality F: 



1)1 



iSr, 
fc 



.0% 



for 77 ^ 1/fc, 



fc-1 



jrff ((l-, 7 )iV-77+-) 



+ 



for r\ > 1/fc, (27) 



with n' = (1 — ??)iV — 77 + 1/fc [see the supplementary 
material for the proof]. The associated bounds for C(S^) 
are obtained by substituting the above expressions into 
Eq. @ together with the identities (|2|). In Fig. |(b) 
we report the one associated to (|27|), together with the 
bounds (|]) and ([m|) , for a direct comparison with Koenig 
and Smith's inequality j| , and with the lower bound (Q) . 
In the low-noise regime of small ./V a capacity bound was 
presented also in Q . 

Bounds for the Amplifying channel A^-' — We now 
prove that the capacity of the channel satisfies in- 
equality (^|) and 

C(A%)^g(KN/[N( K -l) + K }), for TV < ^ , (28) 

see Figs. |l](b) and |](c). The proof of inequality (||) can 
be obtained by adapting the derivation of ( p"3| ) provided 
in Ref. [|l7|. Specifically, the Husimi distribution of a 
generic input state p for m channel uses is 



where |d?) is a to- mode coherent state and o~(a) := 
J ^£ D{p) etf-fl-iF-F-rf-n/i. The corresponding 
output state can then be expressed as [-4^ r ]® m [p] = 
/ d 2m dQ{a) [A%]® m [(x(a)}, while the concavity of the 
output entropy implies 

S{[A N K r m [p])> [ d 2m dQ(d)S([A^ m [a(a)}), (30) 



which is meaningful only when [_4^]® m [o-(d<)] is a quan- 
tum state, i.e. if N(k - 1) ^ 1, when [A%]® m [(x(a)] 
is an TO-mode thermal state with average photon num- 
ber N(k - 1) - 1 per mode, so that S([A^]® m [a(a)]) = 
mg(N(K — 1) — 1). Substituting it into (J30|) we find 



([A%r m )>g(N(K-i) 



1), for N >^(31) 

which, through (|To[) , implies (||). 

Finally, the proof of the bound (|28| ) uses the concate- 
nation^ = _4^o£° withG = N(K-l)+Ka,ndr] = k/G, 
the fact that the capacity is always degraded under chan- 
nel multiplication, and the fact that C(£°) = g(r/N) |J. 

Conclusions: — We have given upper and lower bounds 
for the classical capacity of important active and passive 
bosonic channels, and we have shown that these bounds 
asymptotically coincide (yielding the actual capacity) in 
the regimes of practical interest, i.e. for low transmissiv- 
ity, high thermal noise, or high amplification. 



Q(a) = (a\p\a)/ir m , p = \ d 2rn a Q(a) a(a) , (29) I 10 
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Supplemental Material 

Here we provide explicit derivations of the inequalities ( |2fj| ) and (^7|) . 

Proof of Eq. (p^j: — This inequality can be easily obtained by generalizing to the multimode scenario the beam- 
splitter decomposition of the channel £^ detailed in the Appendix Dl of Rcf. Q [the same decomposition was also 
exploited in Rcf. |?]]]. Consider first the case r\ — 1/k with k an integer [generalization to arbitrary n will given 
later] . The basic idea is to express the transformation induced by the map £y k in terms of a sequence of k — 1 beam- 
splitter interactions that couple the incoming signal mode state p with k independent bosonic thermal baths states 
Pth characterized by the same photon number N. As discussed in Ref. H this can be done in such a way that local 
observers located at each of the k outputs of the array will receive [up to an irrelevant local unitary transformation] 
the same output signal £^ 3 [p\. For instance for k = 3 this can be obtained by setting the transmissivity of the first 
beam splitter equal to r)i = 2/3 and the second one to r\ 2 = 1/2. The same construction clearly can be applied to 
channel [£-f/ fe ]® m of the m-channel use scenario by repeating the decomposition for each channel independently. An 
example of the resulting scheme for k = 3 and m = 2 is shown in Fig. || here A\ and A 2 represent the two channel 
inputs that in principle can be loaded with a non separable state p; B\, C\ are instead the two thermal bath modes 
needed to represent the first channel use, while B 2 and C 2 are those associated with the second channel use [all of 
them being initialized in thermal states having average photon- number N]. In this extended configuration one can 
easily verify that the two-mode states at the ports A\A' 2 , B[B' 2 , and C[C 2 of the figure are all unitarily equivalent to 
the density matrix [f 1 / 3 ]® 2 (/o) [in other words, up to local unitary transformations, each one of those output couples 

yields a unitary dilation || of the same channel [f 1 / 3 ]® 2 ]- This in particular implies that the associated output 
entropies must be identical, i.e. S(A' 1 A' 2 ) = S{B[B' 2 ) = S(C(C^) = S({£^ /3 ]® 2 (pj). Exploiting the sub-additivity of 
the von Neumann entropy fl we can hence write 

S{A' 1 A' 2 B' 1 B' 2 C' 1 C' 2 ) < iS([£^f 2 {p)) , (32) 

where S (A' 1 A 2 B' 1 B 2 C' 1 C 2 ) is the entropy of the joint state at the output of the device. Observing that the transforma- 
tion [i.e., the beam-splitter couplings] that takes the input modes of the system A\A 2 B\B 2 C\C 2 to their associated 
output A[A 2 B[B 2 C' 1 C 2 configuration is unitary, we can then identify S{A\A 2 B' l B' 2 C\C' 2 ) with the input entropy 
S{A\A 2 B\B 2 C\C 2 ). The latter can easily be computed by noticing that the incoming state is just a tensor product of 
p with m(k — 1) = 4 bosonic thermal states with mean photon-number N, i.e., S (AiA 2 BiB 2 C\C 2 ) = S(p) + 4g(N). 
Substituting this into Eq. (^2|) we finally get 

S([^ /3 f 2 (p)) > S(p) + p(N) > ^g(N) . (33) 

The same argument can be easily repeated for arbitrary m and k integers: in this case, we use m(k — 1) local bath 
modes organized in m parallel rows, each containing k — 1 beam-splitter transformations whose transmissivities are 
determined as in Ref. ||. Similarly to the case explicitly discussed above, an inequality for S([£ i/ k ] tg ' rn {p)) can be 
obtained via sub-additivity by grouping the mk output modes into k subsets of m elements each. The resulting 
expression is 

S([£? /k ]® m (p)) > m g(N) . (34) 
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FIG. 3: Beam splitter-decomposition scheme for the channel [£^]® m with 77 = 1/3 and m — 2. Thermal states of mean 
photon- number TV are injected at the input ports Bi, Ci, B2, and C%. 
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FIG. 4: Beam-splitter array used to derive Eq. (J27|) , depicted for m = 2. 



Generalization of this inequality to i] $J 1/fc can finally be obtained along the same lines used in Ref. || by exploiting 
the following composition rules 

KT m °KT m = K't, 2 }® m > ( 35 ) 

which is a trivial multi-mode generalization of the identity (19) from || [here N' = [772(1 — T)i)Ni + (1 — 772)^2]/ (1 — 
771772)]. The reader can check that the resulting expression coincides with the first part of the inequality ((2^). Similarly, 
Eq. (^) can be used to induce a bound for 77 ^ 1/fc by following the same line of reasoning presented in [0 while 
exploiting the composition rule 

[£^f m = [£ r ff m o [Al h ,r m , for 77 ^ rf, (36) 

which is the to- mode counterpart of the identity (B3) of 0. The resulting inequality yields the second part of (f26|). 

Proof of Eq. — The m — 1 version of this inequality was derived in 0, by exploiting a beam-splitter array 
obtained by applying a channel Af n at each of the output ports of the scheme used to derive the m = 1 equivalent 
of Eq. ( |26] ) [see Fig. 12 of [|]. As for Eq. @, the main difficulty in applying the same argument to arbitrary m 
is generalizing such an array to the multi-mode case scenario and properly grouping the corresponding output ports. 
This can be done as sketched in Fig. |j: i.e., adding Af n to each of the ports in Fig. || and by keeping the same grouping 
scheme as before. With this guidance the reader can now closely follow the same derivation given in Ref. (^j [the steps 
are rather cumbersome, but basically coincide with those w e have discussed in the previous section] . 



